Quantum ergodicity on graphs 
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We investigate the equidistribution of the eigenfunctions on quantum graphs in the high-energy 
limit. Our main result is an estimate of the deviations from equidistribution for large well-connected 
graphs. We use an exact field-theoretic expression in terms of a variant of the supersymmetric 
nonlinear cr-model. Our estimate is based on a saddle-point analysis of this expression and leads to 
a criterion for when equidistribution emerges asymptotically in the limit of large graphs. Our theory 
predicts a rate of convergence that is a significant refinement of previous estimates, long-assumed 
to be valid for quantum chaotic systems, agreeing with them in some situations but not all. We 
discuss specific examples for which the theory is tested numerically. 
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Quantum graphs may be viewed as networks of quan- 
tum wires or waveguides and are standard models for 
large molecules and complex quantum systems. Since 
the seminal paper [l[ of Kottos and Smilansky, they have 
become a paradigm for quantum chaos and universal in- 
terference effects which can be observed in the statistics 
of their spectra and of their eigenfunctions (see Q for 
a recent review). Kottos and Smilansky showed (nu- 
merically) that the spectral statistics for a large class 
of quantum graphs follow the predictions of the Gaus- 
sian random-matrix (RM) ensembles up to deviations 
which become smaller when the size of the graph in- 
creases. They also derived an exact trace formula, anal- 
ogous to Gutzwillcr's trace formula which is a scmiclas- 
sical approximation for classically chaotic systems. Sub- 
sequently it was found that some special graphs, namely 
star graphs, exhibit intermediate spectral statistics which 
do not converge to any of the RM expectations 0] . This 
discovery prompted the question as to when the statisti- 
cal properties of the spectrum are describable by means 
of RM theory. This question was answered in |4|]. The 
criterion for RM statistics in the spectrum of quantum 
graphs relies on a well-known quantum-to-classical cor- 
respondence (the corresponding classical dynamics being 
a Markov process on the bonds of the graph) and it is 
generically met in large well-connected graphs. These 
works provide a relatively complete understanding of the 
eigenvalue statistics of quantum graphs. Substantially 
less is known about the eigenfunctions. In this Letter, we 
investigate when the modulus of the high energy eigen- 
functions is uniformly spread (i.e. equidistributes) over 
the graph - a property known as quantum ergodicity. 

It is generally believed, and proved in many cases 
[S IE > tnat an ergodic classical system admits a quan- 
tum ergodic counterpart in the semiclassical limit. De- 
termining the rate at which the quantum ergodic limit is 
approached is one of the major open problems in Quan- 
tum Chaos. General conjectures for this rate have been 
put forward based on heuristic semiclassical arguments, 
RM theory, the random wave model and periodic or- 



bit theory [8[. Surprisingly, establishing when quantum 
graphs are quantum ergodic has proved to be a signifi- 
cant problem. The issue is subtle: Neumann star graphs 
are known not to be quantum ergodic (even though their 
classical dynamics is ergodic) [9(, but quantum ergodic- 
ity has been proved for a family of graphs constructed 
from certain one-dimensional maps [10 ]. Applying the 
field-theoretical approach of Gnutzmann and Altland Q 
enables us to tackle the quantum ergodicity problem for 
graphs, yielding a criterion for when quantum ergodicity 
can be expected for sequences of graphs in the limit as the 
size tends to infinity and, most significantly, an explicit 
estimate of the rate at which the limit is approached. 
Interestingly, in some situations this estimate coincides 
with the expression conjectured in but in others it 
differs from it in a fundamental way. 

In the scattering approach a quantum graph can be 
defined by a pair (G, 5), where G is a metric graph 
with B bonds of finite lengths Lb > (and coordinates 
< Xb < Lb), and S is a 2B x 2B unitary matrix that 
will be discussed below. We assume that the graph is 
simple (at most one bond between any two vertices and 
no loops) and connected, and that the bond lengths are 
incommensurate (rationally independent) and all in some 
fixed interval < L m \ n < Lb < L max . Let us define the 
2B double indices a = (b, d) for directed bonds where 
b = 1, . . . , B refers to the bond and d = ±1 to the two 
directions on each bond (the lower case greek letters a 
and 7 will be used to denote directed bonds through- 
out). The matrix S a . a > contains the probability ampli- 
tudes for a particle to be scattered from a' = (b' , d 1 ) 
to a = (b,d). If the end of a' and beginning of a are 
different vertices then S a ,a' — 0. We assume that all 
scattering amplitudes in S are energy-independent. A 
graph is time-reversal invariant if the scattering matrix 
obeys S = ( j° fl \f ) S T ( ° B ^ ) such that the scatter- 
ing amplitudes for the processes (6, d) — ► (b',d') and 
(6', —d!) — > (6, — d) coincide. The quantum dynamics 
consists in the propagation of one-dimensional waves on 
the bonds of G, with scattering at the vertices encoded 
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by S. The complete quantum dynamics is then encapsu- 
lated in the quantum map defined by the unitary matrix 
U(k) = T(k)S where the diagonal matrix T(k) = e lkL , 
with L = dia,g(Lx, . . . ,Lb,Li, . . . ,Lb), contains the 
phase factors for propagation with wavenumber k > 
from one end of a bond to the other. U(k) describes 
the succession of scattering at a vertex and propagation 
to the next vertex of a wave function ^(x) = {ipb{xb)} 
on the graph, where the wave function is expressed in 
terms of 2B complex amplitudes arb,d) (which we will 
combine into a 2B-dimcnsional vector a) of plane waves 
such that i> b (x b ) = J? (a (f) , + i)e 1 ' !(ll - Lt) + o (( ,,_ 1 )e- 1,:i1 ) 
(Af is a normalization constant). The quantization con- 
dition on a graph implies that the quantum map does 
not change the wave function: U(k)a = a. Equivalently 
U(k) has an eigenvalue unity - a condition that is sat- 
isfied for a discrete set of wave numbers k n which form 
the spectrum. Eigenstates will be denoted by ^ n and the 
corresponding coefficients by the vector a„. 

The classical counterpart of the quantum graph (G, S) 
is a Markov process on the directed bonds of G which 
is obtained by replacing the quantum scattering ampli- 
tudes by the classical probabilities M QjCt / = \U a . a '(k)\ 2 = 
|<Sa. Q '| 2 . The bistochastic matrix M does not depend 
on the wave number and propagates probability dis- 
tributions on the directed bonds [H, 0|. For a con- 
nected graph the classical map M is generically mixing 
(i.e. Um n _ >00 (M") a)0 ,' = ^g-) [11| and always ergodic in 
the stochastic sense. The spectrum of M is confined to 
lie on or within the unit circle with (exactly) one eigen- 
value unity which corresponds to equidistribution on the 
bonds. We write the eigenvalues of M as 1 — with 
mi = and Rcm^ > for i = 2, . . . , 2B - this nota- 
tion anticipates the interpretation of the quantities rrii 
as masses in a field-theoretic setting: massive modes cor- 
respond to decaying modes of the classical map while the 
invariant measure is masslcss. The classical spectral gap 
A g = minj = 2,...2s|wi| determines the slowest decay rate 
in the (time averaged) classical dynamics. It will turn out 
that a simple criterion for quantum ergodicity can be ex- 
pressed in terms of the spectral gap A g , while a more 
detailed criterion uses the full set of masses rrii. Note, 
that the topology of a graph influences the spectrum, 
e.g. any graph which can be cut into two disconnected 
components of similar size by just erasing one bond has 
a small gap which is at most of order 1/B. 

Our concern here is to determine how the prob- 
ability density associated with the eigenstate ^„ is 
spread over G as k„ — > oo. Let us introduce observ- 
ables that are constant on each bond. These can be 
represented by real diagonal 2B x 2B matrices V = 
diag(Vi, . . . , Vb, Vi, . . . , Vg). The expectation value of 
the observable V in an eigenstate *ff n is given by 

<y) n =5> f Lb d Xb \u^)\ 2 = ^"i^i a ») + 0( fc -i) . 



With the spectral counting function N(k) — J2 kn 9{k — 
k„) ~ ^2^-k the spectral average of the expectation value 
of the observable V is given by the metric average over 
the graph 

v ; k n <K 

This statement is known as the local Weyl law and can be 
recovered for any finite graphs straightforwardly, e.g. by 
the periodic-orbit approach developed in For a uni- 
formly distributed eigenstate the coefficients a a are all 
of equal modulus such that the expectation value {V) n 
coincides with the spectral mean for any observable V. 
Quantum unique ergodicity is the statement that such 
uniform distributions are obtained in the high energy 
limit of any subsequence of states \E' rij with n,-_|_i > n,, 
i.e. one has Um,-_ >00 (V) ni = Ay. This property turns out 
to be too strong for graphs and cannot be expected to be 
realized because it is known that sequences of states exist 
that are scarred by short periodic orbits 0, [12] ■ The less 
strong property of quantum ergodicity can be defined by 
the vanishing of the variance 

i.e. Fy = for all observables V. This implies the ex- 
istence of subsequences of eigenstates of density one for 
which limi^ 00 (V) ni — Ay. It will be seen later, see ([6]), 
that Fy in general does not vanish for a finite graph. In- 
deed one may expect, as for spectral statistics, that quan- 
tum ergodicity can only be attained in the limit B — > oo. 
One therefore considers a sequence {(G;, Si)} of quantum 
graphs having increasing numbers of bonds -B;+i > 
For each I, one defines Aiy t and Fiy t by formulas {J) and 
([2|) applied to (G;, Si). One still has to specify the accept- 
able sequences {V;} of observables. For these it is natural 
to require the existence of lim^oo Aiy l and a bound from 
above for |V; jQ ,|. Then a sequence {(G;, Si)} is said to be 
asymptotically quantum ergodic if -Fz.y, — > as I — > oo 
for any acceptable sequence of observables {V;}. In order 
to avoid cumbersome notation, we will separately work 
on a single quantum graph (T;, S/) of the sequence and 
drop the index I. Without loss of generality we will also 
assume A\_y t = 0. 

Our aim is to show that a large class of such sequences 
exhibit asymptotic quantum ergodicity and to give ex- 
plicit conditions which distinguish between ergodic and 
non-ergodic sequences. In the course of our derivation we 
will also obtain an estimate of Fy for large (but finite) 
quantum graphs. While previous attempts have largely 
relied on periodic-orbit theory we will here employ a field- 
theoretic approach. In order to implement our approach 
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let us first introduce 



Fy = lim — — V (V)l 



k n <K 
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a,a' — 1 



which vanishes asymptotically if and only if i*V vanishes 
(note that we assume Ay = 0). For all practical pur- 
poses, Fy and Fy are the same for a large graph. In ((3|) 
we have implicitly defined £ QQ ,< which will be the central 
quantity of interest. It can be expressed as 
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(4) 



dct(C/ e (fc)t 



where (. . .) k = lim^oo J Q 



K 



dk, U e (k) = 1 - 

e~ e U(k) and = <Va"W " 

Using a standard procedure, £ Qa < can be expressed ex- 
actly in terms of a variant of the supcrsymmetric non- 
linear a- model Q ■ The procedure has two main ingredi- 
ents: (i) for incommensurate bond lengths the integral 
over k in ^ can be replaced exactly by an integral over 
B independent phases fl3j |. (ii) Zirnbauer's color-flavor 
transformation [141 ]. For broken time-reversal symmetry 
the (7-model for the graph leads to the exact expression 



hm e — — — 

e ^o+ dj + dj_ 
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(5) 



d(Z, Z) sdct(l - ZZ) sdet _1 (l - ZS+ZS- 



Here Z = diag(Zi, . . . , Z B ) and Z = diag(Zi, . . . , Z B ) 
are 4B x 4B block-diagonal supermatrices where each 
block is a 4 x 4 matrix Zt, t dd',ss' (d — ±1 is the direc- 
tion index and s = B, F refers to bosonic and fermionic 
sectors of the supermatrix, see 0] for further details). 
The matrices S+ = (1 - j + E a ' a P FF )S and Si = 

5^(1 — j-E aa Pff) (where Pff is the projector onto 
the fermionic sector) contain the scattering amplitudes 
and source terms oc j±. For time-reversal symmetric 
graphs an exact cr-model can be derived with slightly 
more technical effort. The result looks formally simi- 
lar to (|5|) with the superdetcrminants replaced by their 
square roots and with the dimension of the matrices dou- 
bled (see 0]). Here we will give explicit formulae only 
for broken time-reversal symmetry but will state the fi- 
nal results also for the time-reversal invariant case. More 
details will be given elsewhere . 

We calculate the integral over Z and Z in a 
saddle-point approach (writing sdet ±1 (l — ZAZB) = 



cxp(±str log(l — ZAZB))) where the saddle-point anal- 
ysis is performed at j± = 0. The saddle-point equa- 
tions lead to a mean-field solution Z&^'.ss' = Y ss >5dd', 
ZbAd'.ss' = Y ss '5dd' where the 2x2 matrices Y and Y 
span the saddle-point manifold. Note that the depen- 
dence of the integrand in ([5]) on the scattering matrix 
S of the graph drops out for mean-field configurations so 
that these give a universal contribution to £ aa i . Configu- 
rations orthogonal to the saddle-point manifold are taken 
into account in a Gaussian approximation - here, cru- 
cially, the system dependence does not drop out. There is 
a direct correspondence between the modes of the classi- 
cal map M and the field configurations. Equidistribution 
on the graph corresponds to the unique eigenvalue unity 
of M on one side and to the mean-field configuration on 
the other side. This is the masslcss mode which is calcu- 
lated exactly in our approach. The decaying modes of the 
classical map correspond to the (massive) configurations 
of the supermatrices Z and Z which are orthogonal to the 
saddle-point manifold. The contribution of the massive 
modes (and their coupling to the mean-field mode) are 
only taken into account approximately. This goes beyond 
the diagonal approximation in a periodic-orbit approach 
to £ QQ ' , which is equivalent to a Gaussian approximation 
of all modes in the field-theoretic approach and which ac- 
tually diverges like 1/e for the quantity we are interested 
in (further details will be published elsewhere [IH). 

The contributions from the mean-field configurations 
and the Gaussian fluctuations turn out to be additive 
EH and thus we may write the fluctuations as Fy = 
Jy auss . An explicit calculation yields 
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for the universal contribution. Here j3 = 1 for broken 
time-reversal invariance and (3 = 2 for time-reversal in- 
variant graphs. The integration over the Gaussian fluc- 
tuation yields 



rriGauss 



(trL) : 
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2ta2 



1 — M 



-L'V 



(7) 



where the trace is over all massive modes such that the 
unit eigenvalue of M is excluded. 

Equations © and Q arc our main results. Crucially, 
([7]) shows that the behavior of Fy for large graphs de- 
pends on the spectrum of the classical map M. Slowly 
decaying classical modes with eigenvalues near 1 give rise 
to large deviations from universality in the rate of quan- 
tum ergodicity. This does not conform to expectations 
based on the conjectures in 0], which correspond to ig- 
noring the Gaussian fluctuations. 

In the remainder we discuss the different behaviors 
that sequences of quantum graph may exhibit depend- 
ing on the spectral properties of M and on the range 
of validity of equations ^ and ([7]). The dependence of 
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FIG. 1: (color online) Deviations Fy from quantum ergod- 
icity for sequences of complete graphs. Data points refer to 
numerical calculations and lines refer to the saddle-point ap- 
proximation. From top to bottom: complete Neumann graph 
(/3 = 2. red), complete Neumann graph with magnetic field 
{(3=1, light blue), complete DFT graph (/3 = 2, dark blue), 
complete DFT graph with magnetic field (/? = 1, green). 



Fy on the observable V is of minor importance here and 
we restrict ourselves to the most simple generic choice, 
VbLf, = ±const 



contributions to Fy which may be calculated by going 
beyond the Gaussian approximation. Thus while we be- 
lieve that our criterion for quantum ergodicity A g £> — > oo 
is correct we do not expect the Gaussian approximation 
to give the exact asymptotic result if Fy is dominated by 
the massive contributions. We have checked numerically 
the validity of the saddle-point approximation for two se- 
quences of complete graphs with different scattering ma- 
trices; see Fig. 1. In the first sequence we have used dis- 
crete Fourier transform (DFT) matrices at each vertex. 
In order to break time-reversal invariance we also added 
magnetic fields. In this case l'™ a,liidd > p^ uss w 
and the numerics show perfect agreement. In the sec- 
ond sequence we have used a scattering matrix which 
is equivalent to the so-called Neumann (or Kirchhoff) 
boundary conditions at each vertex. These strongly en- 
hance backscattering which leads to a large number of 
small masses. As one may expect, the saddle-point ap- 
proximation is far from perfect in this case but still gives 
a reasonable estimate of the rate. 

We would like to thank Uzy Smilansky for invaluable 
discussions. 
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Obviously the universal part vanishes in the limit B — > oo 
without any conditions on the graph. We conjecture that 
the sequence of graphs is asymptotically quantum crgodic 
if and only if Fy auss also vanishes in this limit. If the 
spectral gap remains finite as B — > oo the contribution 
from the massive modes remains at most of same or- 
der as the mean-field contribution. If the spectral gap 
decreases as B — > oo we may assume A 9 ~ B~ a with 
a > 0. The contribution from the Gaussian fluctuations 



can then be bounded by 



^Gauss 



< -B- 

- BA„ 



B a ~ x . Based 



on the saddle-point approximation we can thus give the 
following criteria: (i) a sequence of quantum graphs is 
asymptotically quantum ergodic if BA g — > oo i.e. a < I, 
[ii) a sequence of quantum graphs is not asymptotically 
quantum crgodic if a > 2, (ra) if 1 < a < 2 asymptotic 
quantum ergodicity depends on all masses (i.e. on the 
complete spectrum of the classical map). 

Finally, let us comment on the validity of the saddle- 
point approximation in this context. The masses rrii de- 
termine the stability of the saddle-point manifold - for a 
single mode the Gaussian approximation can be expected 
to capture the dominant contribution if rrii ^> 1/B. 
When this condition is not satisfied the saddle-point ap- 
proximation is not valid. However, even if rrii 3> 1/B 
the number of modes still increases proportionally to B, 
which makes it hard to estimate the error. Especially if 
there are many small masses one may expect additional 
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